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Abstract 

We consider random non-hermitean matrices in the large N limit. The 
power of analytic function theory cannot be brought to bear directly to analyze 
non-hermitean random matrices, in contrast to hermitean random matrices. 
To overcome this difficulty, we show that associated to each ensemble of non- 
hermitean matrices there is an auxiliary ensemble of random hermitean matrices 
which can be analyzed by the usual methods. We then extract the Green's 
function and the density of eigenvalues of the non-hermitean ensemble from 
those of the auxiliary ensemble. We apply this "method of hermitization" to 
several examples, and discuss a number of related issues. 



'I found it difficult,...., to keep my mind from wandering 
into the magical world of random matrices. " 
F. J. Dysonft 

1 Introduction 

When Wigner ]T[ first proposed treating the Hamiltonians of complex systems as ran- 
dom matrices, he quite naturally took the matrices to be hermitean (or real symmetric 
for time reversal invariant systems) in accordance with the fundamental postulates of 
quantum mechanics. In the last few years, there has been increasing interest in con- 
sidering non-hermitean random matrices in the context of various physical problems. 
(We will use the term "non-hermitean" to include "real but not symmetric") We 
will mention a few examples: QCD at finite chemical potential^], nuclear decay ||, 
dissipative systems || |J , and neural networks || . In all these examples, the random 
non-hermitean matrices were taken from a Gaussian distribution. Recently, Hatano 
and Nelson || (see also 0) have used non-hermitean random matrices to study the 
pinning of magnetic flux lines in high temperature superconductors. 

In this paper, we discuss some aspects of random non-hermitean matrix theory, 
hoping to clarify some of the existing (and growing, see e.g. [p] for some recent exam- 
ples) literature, which in the authors's opinion, contains a number of confusing and 
confused statements. We develop a method to deal with non-hermitean matrices, 
by reducing the non-hermitean problem to an auxiliary hermitean problem. Unlike 
many examples that appear in the literature, our formalism is generic and is not 
limited to Gaussian ensembles. We then present some specific results. Some of our 
results may already be known in the literature. In particular, as we were completing 
this work, a series of very interesting papers by Zahed and collaborators 0] studying 
non-hermitean matrices appeared. 

A basic tool in studying hermitean random matrices <ft (henceforth all matrices 

1 F.J. Dyson, "Selected Papers of Freeman Dyson with Commentary" 
(American Mathematical Society, International Press, 1996), p. 39. 
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will be taken to be N x iV with N tending to infinity unless otherwise specified) is 
the Green's function defined by 

G < 2 > = <>I^> M 

where < • • • > denotes averaging over the random distribution from which the ma- 
trices <fi are drawn. Diagonalizing by a unitary transformation we have 

where the N real numbers {A^} are the eigenvalues of <p. Thus, G(z) is a meromorphic 
function with poles at the eigenvalues of </>. In the large N limit, the poles merge into 
a cut (or several cuts) on the real axis. Powerful theorems from the theory of analytic 
functions can then be brought to bear to the problem of determining G(z) |10| . All 
of this is well-known. In contrast, for a non-hermitean matrix, the eigenvalues invade 
the complex plane. For example, for non-hermitean matrices <fr generated with the 
probability P{4>) = 4e _iVtr< ^, Ginibre determined long ago JO]] that the eigenvalues 
are uniformly distributed over a disk of radius unity in the complex plane. Another 
simple example is presented in the Appendix. Thus, we lose the powerful aid of 
analytic function theory. To get oriented, consider the simple example of the function 

d9 1 f dw 

9{z) ' 



2tt z — e %e J 2iri w(z — w) 

which appears to depend only on z. In fact, g(z, z*) = 9 (\z\ 2 — 1) jz is clearly non- 
analytic, and depends on both z and z*. Speaking mathematically, one would say 
that this integral defines two functions, one for \z\ > 1 and one for \z\ < 1. 

In the theory of hermitean random matrices, a powerful method consists of the 
Feynman diagrammatic expansion ||12|| in which one expands G(z) = Z)fc(tr0 fc ) / z k+1 
as a series in 1/z, the " bare quark propagator." This method is clearly no longer 
available in studying non-hermitean random matrices. The knowledge of G(z) as a 
series in 1/z can no longer tell us anything about the behavior G(z) for small z, as 
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the simple example above makes clear. The eigenvalues fill a two-dimensional region 
rather than a one-dimensional region, as is the case for hermitean matrices. 

A renormalization group inspired method used in T^, 13, W, |TIJ implicitly involves 
an expansion in 1/z and is thus also not available for dealing with non- hermitean 
matrices without suitable further developments. 

This paper is organized as follows. In Section 2 we review some of the basic formal- 
ism. We also make a tentative connection with the theory of quaternions. Section 3 is 
the heart of this paper which consists of showing that the Green's functions and den- 
sity of eigenvalues of the random non- hermitean matrix may be determined from those 
of an auxiliary hermitean matrix. The hermitean auxiliary ensemble may be analyzed 
by the usual methods. We refer to this process as the "Method of Hermitization" . 
In Section 4 we show how to determine the boundary of the eigenvalue distribution 
of the complex matrix without ever determining its density of eigenvalues p(x,y). In 
Section 5 we present some simple applications of the hermitization method. In Section 
6 we discuss some aspects of "energy-level" dynamics of non-hermitean hamiltonians. 
In particular, we show how the celebrated phenomenon of level repulsion in the case 
hermitean hamiltonians changes when one discusses non-hermitean hamiltonians. We 
present a simple (yet instructive) example in the Appendix. 
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2 Some Basic Formalism 



In this section, we will derive the central formula for studying the eigenvalue dis- 
tribution of random non-hermitean matrices. While this formula is already known 
in the literature, we hope to clarify a number of misleading statements made in its 
derivation. In order to be absolutely clear, we will begin with some almost laughably 
elementary facts. 

For z = x + iy we define 

o d 1 f d . d\ 

d =d- Z = 2{d X --%) (2 ' 1} 
so that dz = 1. Clearly, dz* = 0. (We denote complex conjugation by *.) It follows 
that dz* n = (n > 0), and d f(z*) = in a region in which / has a series expansion 
in powers of z* . As is well-known, this fails if / does not have a series expansion. In 
the simplest example, we have 

d L = nS(x)S(y) , (2.2) 

z 

which we can easily prove by integrating the left and right hand sides over a small 
square centered at the origin. Similarly, we define 

so that d* z* = 1. We have, obviously, d* (1/z) = ir5(x)5(y). 

Let us now calculate dd* log (zz*) = dd* (log z + log z*). At this point it is 
important to keep in mind that d log z is not \ in the cut plane! [d log z is \ on the 
multi-sheeted Riemann surface of course.) Instead, the correct formula is 

d log z = - + 7i6(-x)5{y) (2.4) 

z 

with the standard choice of cutting log z along the negative real axis. Thus, 

dd* (log z + log z*) = d*(- + ir9(-x)6(y)) + c.c. 

z 



= ir5(x)5(y) + n5(x)5(y) - n5(x)5(y) , (2.5) 
4 



where the last term comes from (d + d*) ir6(—x)5(y) = —7r5(x)5(y). We finally obtain 



dd* log (zz*) = 7r5(x)S(y) . (2.6) 

Notice that if we had glibly used d log z = -, as was sometimes done in the literature, 
we would have missed the last term in (|2.5|) and obtained an erroneous formula off by 
a factor of 2 from the correct formula, thus leading to an endless stream of apparent 
paradoxes. 

With these basic facts established we are now ready to deal with the average 
density of eigenvalues of a non-hermitean matrix <fi. We diagonalize the matrix by a 
similarity transformation 

<P = S^AS (2.7) 

where A denotes a diagonal matrix with elements Aj, % = 1,...,N. Taking the her- 
mitean conjugate we have <ft = S^A*S~^. By definition, the density of eigenvalues 
is 

p(*> y) = (jf E 5 ( x - Re A *) <% - Im A *)> • ( 2 - 8 ) 



Given the preceding discussion, we obtain from (O) that d* tr {z — A) 1 = d* tr (z 



0) 1 = ?r J2i $( x ~ R- e K) d{y ~ I m \)- Thus, using the identity 

det (z - A) (2* - A*) = det {z - <f>)(z* - <fl) 
we may express p in terms of as 

p(x,y) = -dd* (itr log (* - 0)(z* - 0t)) . (2.9) 

7T iV 

This formula has been known in the literature (see e.g., [H]), but there a seemingly 
unnatural small regulator had been introduced to make the argument of the logarithm 
in ( p.9|) strictly positive to avoid the branch point at 0. Our derivation of ( p.6|) suggests 
that that regulator is not really necessary. 



Our expression (2^9) is manifestly symmetric in z and z*. It arises quite naturally 



in calculations involving Grassmann variables^, [| where det (z — (f))(z* — <ft) is the 
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"fermion determinant" . However, it involves the logarithm inside the average, as well 
as two derivatives. Thus, unless one has a simple way of calculating the average in 
( |2.9| ) (e.g., the replica method in the case of Gaussian ensembles ||), ( |2.9| ) is not 
a practical expression for p(x,y). In most cases it is easier to calculate the Green's 
function 

G{z,z*) = (^tr— *— > = (i V- — ^- ->, (2.10) 

(with Xi = Xi + iyt) than to do the average (^tr log (z — 4>){z* — $)) . Then, from 
( |2.8| ) and from the complex-conjugate of ( |2.2| ), we have 

p(x,y)^-d* G(z,z*), (2.11) 

7T 



which is a simpler representation of p(x, y) than (J2T 

These two representations for p(x, y) and the relation between them can be in- 



terpreted by recognizing (|2.9|) as a two dimensional Poisson equation for the elec- 
trostatic potential — ~(4 tr log [z — 4>)(z* — $)) created by the charge density 
p(x,y). This connection between eigenvalue distributions of complex matrices and 
two-dimensional electrostatics has been long known in the literature || ||. Continu- 
ing along this line, consider the Green's function G(z,z*) in ( |2.10| ). If we define the 
electric field E = (Re G, — Im G) then from the definition of G(z, z*) we have 



E{x) = I d 2 y p(y)T^- 2 (2.12) 

f ~ y\ 



and thus ( ^i.ll ) is simply the statement of Gauss' law for this electrostatic problem. 



We end this section with a speculative remark. As we pointed out earlier, in the 
theory of random hermitean matrices it is well known that it is much easier to work 
with G(z) = (jf^-^z) rather than with the density of eigenvalues p(p) directly, since 
the power of analytic function can be brought to bear on G(z). Of course, one can 
go from p(x) to G(z) with the identity 

-lm^— = -^—^5(x). (2.13) 

7T X — It 7T x z + e z 

Here, for random non-hermitean matrices our formula ( |2.9|) gives the density of 
eigenvalues p(x, y) directly, but this formula is awkward to work with because of the 



logarithm. Its direct analog in the theory of random hermitean matrices would be 
p(x) = -^{jjti6{x — 0)), which would also be extremely awkward to work with. It 
would thus be desirable to define a quantity analogous to G(z) and develop a method 
for calculating it. Reasoning along these lines we are led to an attempt to write a 
function of a quarternionic variable, in the same way that going from p(x) to G(z) we 
went from a function of a real variable to a function of a complex variable. Indeed, 
the obvious analog of Q2.13Q is easy to find, namely, 



1 

2j 



{z-je)\z- je\ 1 (z - je)\z - je\ 



je\ 3 ( X 2 + y 2 + £ 2)3/2 



2nS{x) 6(y) . (2.14) 



Here z = x + iy is a complex number, e is a small positive real number and {1, i, j, k} 
is the standard basis of the quaternion algebra^]. (The absolute value of a quarternion 
is defined by \a + ib + jc + kd\ = y/a 2 + b 2 + c 2 + d 2 .) 

Unfortunately, ( p,14| ) as it stands, is less useful than (|2.9| ), because it leads to the 
quaternionic Green's function G(q) (q being a quaternionic variable) 

which involves the absolute value operation explicitly, and thus cannot be written 
as a simple trace like ( |1.1| ). Note, however, that if we split q quite generally into 
q = z + jw (z,w being ordinary complex variables), then G(z + jw) is manifestly 
non-analytic in z as w — > 0, even before taking the average. 

We pose this as an interesting problem, namely to find a useful quaternionic gen- 
eralization of ( |1.1|) for random non-hermitean matrices, and to develop the analogue 



of the work of Brezin et al. 10 associated with it. 



2 One may of course write an equivalent formula with j replaced everywhere by the third quater- 
nion basis element k. Also, e may be taken complex. 
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3 The Hermitization Method: Reduction to Ran- 
dom Hermitean Matrices 



We have emphasized that a straightforward diagrammatic method is not allowed for 
non-hermitean matrices. Somewhat remarkably, we can arrive at a diagrammatic 
method indirectly. We will show that the problem of determining the eigenvalue den- 
sity of random non-hermitean matrices can be reduced to the problem of determining 
the eigenvalue density of random hermitean matrices, for which the diagrammatic 
method may be applied. 

We start with the representation p(x, y) = ^dd* (-^trlog (z— 6)(z*— <fft)) (namely, 
Eq. ( |2.9D . ) By standard manipulations we observe that (tr log (z — 6)(z* — $)) = 
(tr (2JV) log H) — inN 2 , where H is the hermitean 2N x 2N matrix Hamiltonian 

6-i 



H = 

Thus, if we can determine 

F( m z,z*) 



(3.1) 



2N 







( tr ( 2iv) lo S (V-H)) 



(3.2) 



we can determine p(x,y). 

Consider now the propagator associated with H, namely, 

1 



GStoz,**) = (I 



J/ 



H 



) 



(3.3) 



where rj is a complex variable and the indices p and v run over all possible 2iV 
values. Here we followed a common practice and borrowed some terminology from 
gauge field theory: we may consider 0, 6' as "gluons" (in zero space-time dimensions), 
which interact with a 2N dimensional multiplet of "quarks" ip 11 , with a complex mass 
matrix (the "inverse propagator") 

( r\ z\ 



Go 1 



V z* 77 J 



(3.4) 



(expressed in terms of its N x N blocks.) The crucial point is that since H is her- 
mitean, Q£ can be determined by the usual methods of hermitean random matrix 
theory. In particular, as we already mentioned, the diagrammatic evaluation of ( |3.8| ) 
is essentially the expansion of Q£ in powers of I/77, with interaction vertices H. This 
is a well defined procedure for large 77, and it converges to a function which is analytic 
in the complex rj plane, except for the cut (or cuts) along the real axis which contain 
the eigenvalues of H. After summing this series (and thus determining Q^{f]] z, z*) in 
closed form), we are allowed to set r\ — > in (|3.3|) . Speaking colloquially, we may say 
that the crucial maneuver here is that while we cannot expand in powers of z, we can 
arrange for z to "hitch a ride" with r\ and at the end of the ride, throw 77 away. 

We can now calculate p(x, y) by two different methods, each with its advantages. 
The first is to simply observe that 



1 

H 



( 



V A 

Thus, the quantity is simply the lower left block of G£(rj = 0;z,z*). In other 
words, once we have Q£(r);z,z*) we can set 77 to zero and use ( j2.10| ) and ( 2.11 ) to 



write 

' 1 N 




(3.5) 



This observation is the basis for many of the calculations in 0. We will illustrate 
this discussion with a simple example below. 
An alternative is to take the trace of 

^^*) = ^V< tr c-> ^> = (3 - 6) 



from which one can determine 

We refer to this procedure of calculating G(z,z*) = ^(tr^r) and p(x,y) by 
manipulating the hermitean matrix H as the "Method of Hermit iz at ion" . 

We now derive a dispersive representation of F(r];z,z*) in terms of G(i]; z, z*). 
Recall that the eigenvalue density of H is given by the discontinuity of Q accross the 
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cut, namely. 

11 

u(n',z,z*) = — (tr (2JV) S(n-H)) = -lm.g{fi-ie)\ e ^ 0+ . (3.7) 

From (|2.4j) we then find 

^-Ffa; z, z*) = ^(77; -z, z*) + vr [1 - fi(Re 77)] 5(Im 77) , (3.8) 

where SI (Re 77; z, z*) is the integrated density of eigenvalues of H, 

S7(/W*) = tr (2JV) Ofr-H)). (3.9) 

As a result of the "chiral" block structure of if, its spectrum consists of pairs of 
eigenvalues with equal magnitudes and opposite signs. Namely, if £ is an eigenvalue 
of H, so is — £. This symmetry of the spectrum of if around /x = Re = implies 
S7(0) = 2- Thus, setting rj = is pure imaginary we obtain 

— z, z*) = iQ{is- z, z*) + -5(s) . (3.10) 
as 2 

Integrating over s we have 

s 

F(i0+; z, z*) = F(is; z, z*) - i J dsg(is; z, z*) (3.11) 

0+ 



where the lower integration limit avoids the 5(s) piece in ( p. 10 ). From the definition 



of F, we have that F(rj; z, z*) — > log 77, independent of z, as 77 tends to infinity. Thus, 
applying the operator |<9<9* to both sides of ( |3.11| ), we finally obtain 

00 

p( X) y ) = / dsdd*g{is; z, z*) . (3.12) 

7T J 

0+ 

Thus, the density of eigenvalues of the non-hermitean random matrices <p is deter- 
mined by an integral over Q(r];z,z*). The determination of g(i];z,z*) involves only 
the well-known and more elementary problem of determining the density of eigen- 
values of the hermitean random matrices H, which is the essense of our method of 
hermitization. 
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In fact, H in 



ID consists of the sum of of the deterministic matrix 
-z\ 



Ho 



J 



(3.13) 



and the random matrix 



V 



( ° 

V 4 



(3.14) 







Note that the deterministic piece has N eigenvalues +r = \z\ and N eigenvalues 
— r. This problem of adding a deterministic matrix and a random matrix has been 



extensively discussed in the literature [12], [T7|, |T|, H3 |20|, The Green's function 
( |3.6| ) and the density of eigenvalues ( |3.7| ) are related by the standard expression 

no OO 

u(n;z,z ) f . ^ 2r/ 



£(77; 2, 2*) = J 



dfx- 



7}- (i 



d\x u(n\ z, z* 



T] 2 — fl 2 



(3.15) 



Here we used the fact that u(fi; z, z*) = u(—fj,;z,z*), namely, the symmetry of the 
spectrum of H mentioned above. 

Using the "dispersive representation" (|3.15 ) we can derive a compact formal ex- 
pression for the density of eigenvalues of the non-hermitean random matrix 0. Insert- 
ing ( [3.15D into (|3.12|) we obtain 



dfidd*uj(fi; z, z* 



ds- 



s 2 + jJ, 2 



(3.16) 



0+ 

Note that the "+" can now be safely removed from the lower limit of the s integral. 

Also, note that while the s integral appears to be logarithmithically divergent, it 

00 

is in fact convergent thanks to the "sum rule" / d//o;(/i; z, z*) = 1. To make this 

—00 

last observation explicit, let us recall that cu(fx;z,z*) = dn ^^ z ) ; where Q(/j,,z,z*) 
is the integrated density of eigenvalues defined in ( |3.9| ). We can thus integrate over 
\i by parts in ( p.!6|) , thereby bringing in z,z*). We obtain finally the compact 
expression 

n(fi;z,z*) 



dfidd* 



(3.17) 
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To summarize, we have obtained a formalism (the "Method of Hermitization" ) 
for reducing the problem of dealing with random non-hermitean matrices to the 
well-studied problem of dealing with random hermitean matrices. Given the non- 
hermitean matrix 0, we study the hermitean matrix H instead. By whatever method 
one prefers, once one has determined the quark propagator Q$ (or its trace, the 
Green's function Q{j]\ z, z*)), one can in principle obtain p(x, y) using (|Q| ) or (|3.17|) . 
Whether that can be done in practice is of course another story. 

Note that a "shadow" of the variable rj appears in Section 3 of the first paper in||, 
as a small real regulator A. This regulator explicitly "breaks" what these authors refer 
to as the "holomorphic symmetry", by coupling "quarks" and "conjugate quarks", in 
an anlogous manner to the way a small explicit quark mass breaks the chiral symmetry 
in QCD. Then N tends to infinity first, followed by A — > 0. After these limits were 
taken, in that particular order, one observes that (tr(z — 0) _1 ) is not a holomorphic 
function. Since by now A = 0, these authors refer to this phenomenon as "spontaneous 
breakdown of holomorphic symmetry" . From the point of the hermitization method, 
we see that it is rather misleading to treat rj = A as a small regulator. On the contrary 
- one expands in powers of 1/rj for 77 large, and only after summing this perturbative 
series for the hermitean matrix H, one is allowed to set 77 = 0. From this point of view, 
as well as from the simple example given in the Introduction, the non-holomorphy of 
(tr(z — is not more mysterious than the fact that has complex eigenvalues. 

It is worth emphasizing that the formalism developed here has nothing to do with 
large N as such. The formalism is of course indifferent to the method one may choose 
to use to determine Q^,{r\\ z,z*). The problem of determining the Green's function of 
chiral matrices such as H has been discussed by numerous authors [JIB], [2(], 
In some cases, results can be obtained for finite N |22| . 



It is also clear that the formalism here is totally independent of the form of the 
probability distribution P((f),(f)>). In special cases one can say more. For example, 
suppose that P(4>, ft) is invariant under 



e ia (p, ft-^e~ ia ft. (3.18) 
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For instance, this holds if P is constructed out of tr (<^</>) m , (tr <p^<f)) n , etc. Then, 
clearly, (^tr log (z* — (fr)(z — 0)) is a function of \z\ 2 = r 2 rather than of z and z* 
separately. The formalism just given simplifies somewhat: the propagator Q^{j]\ z, z*) 
and the functions Q(r];z,z*), u(p;z,z*), and fl(p; z, z*) do not depend on z and z* 
separately, but only on r = \z\. Acting on such functions, the operator dd* in ( |3.17| ) 



reduces to 

n ^d 9 dldd , 
QQ* — r 2 = r (3.19) 

dr 2 dr 2 4r dr dr 
The density of eigenvalues p(x, y) = p(r) is obviously circularly symmetric. In par- 
ticular, the Gaussian problem considered by Ginibre is of this class. 

In some problems, one might be interested in the density of zero eigenvalues. 
Expanding Q(p;r) = fi (/u) + r2 ^i(A t ) + ••• as a series in r 2 , and inserting into 
(|3.17|) we obtain formally a compact expression for the density of zero eigenvalues 

p(o) = /fn 1 (A*) 

It is worth emphasizing again the point we made in the Introduction, namely 
that we can no longer expand in powers of 1/z and use the diagrammatic expansion 
directly. Consider in particular the class of distribution we just mentioned. Suppose 
we are allowed to expand in powers of -, then we would compute 

p(x,y) = ^dd* (-Ur log (z - 4>)) + c.c. 

= 5{x)5{y) + -dd* i £ ^P- + c.c. = 5(x)8(y) , 
7i Jy £ -r z K 

since in each term (tr (f) k ) vanishes (for k > ) by the symmetry of the distribution 
P(4>, (ft). We would have concluded erroneously that p(x,y) is concentrated in a 
5-function spike at the origin. 

In contrast, it is of course perfectly legitimate to expand in powers of I/77 and 
use the diagrammatic method to compute G£(t]; z, z*) associated with the hermitean 
matrix problem. 

A similar point worth re-emphasizing is that the formula ( 2.1 1| ) for p(x, y) is also 



not very useful as it stands unless we have a way of computing G(z,z*). What we 
have to do is to compute Q^(tj; z, z*) by a diagrammatic method and then use (|3.5| ). 
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Let us illustrate some of this discussion with a simple example, the Gaussian 
problem first worked out by Ginibre|ll|, namely the density of eigenvalues of non- 



hermitean matrices randomly generated according to the normalized probability dis- 
tribution 

P(0) = }_ e ~Ntr4>H _ (3>20) 
Zj 

According to our method of hermitization, we can reduce this problem to that of 
obtaining the density of eigenvalues of the hermitean Hamiltonian H in ( |3.1| ), con- 
sisting of a deterministic matrix and a hermitean Gaussian random matrix of a chiral 
form. As remarked earlier, this belongs to a class of problems long discussed in the 



literature |T|, [17|, [18], [IJ, gg, p. 

Let us first apply the by-now well-known diagrammatic method to calculate the 
Green's function Q(rj;z,z*) = Q(rj;r) in the large N limit. Recall that the quark 
self-energy matrix E£ and the quark propagator Q% in (|3.3|) are related by 

' i V 



(3.21) 



Since S is the sum over all one-quark irreducible graphs we have the Dyson- Schwinger 
identity 

Vt = {VGV)$, (3.22) 
where V was defined in ( |3.14 ). After a simple calculation we find from ( ^.20[ ) and 



(|3722|) that 

K = Q{mr)Si;. (3.23) 



Thus, from ( |3.21| ) and ( |3.23| ) we observe that the diagonal blocks of Q£ are equal 



and proportional to the N dimensional unit matrix. By "index democracy" this 
proportionality constant is nothing but the Green's function G(i];r) = Q. We thus 
have 

' ' ' (3.24) 



F = — 

v r 2 - ( v - gy 



Z* Q Tj 

Tracing both sides of ( |3.24| ) we obtain the "gap equation" 



G 3 -27]G 2 + (l + 7] 2 -r 2 )g -77 = 0. (3.25) 
14 



The relevant solution can be immediately written down using the classic Cardano 
formula. (Of the three roots, the relevant root is the one which goes as I/77 as 77 tends 
to infinity, as is evident from the definition of Q(rj;r).) We will not bother to write 
the expression explicitly Then, following (|3.5|) we find 



p{x t y) = -Sr Z * (3.26) 
Ti r l — fer (0; r) 

Solving ( ggg ) at 77 = we find that either 0(0; r) = or r 2 - £ 2 (0;r) = 1. 
Plugging Q = into (|3.26j ) we obtain an expression proportional to S(x)5(y) which is 
nothing but the zero density outside the eigenvalue distribution. In the other case, 
substituting r 2 — Q 2 (0;r) = 1 in (|3.26|) we obtain Tip = d*z* = 1, which corresponds 
to a uniform p inside the disk of eigenvalues. These two possibilities match at r = 1 
which must therefore correspond to the boundary of the eigenvalue distribution. Thus, 
finally 

'1 , r < 1 

TPotanJfotfJH ( 3 - 27 ) 
, r > 1 



in agreement with JO], but the derivation here is considerably simpler. It is worth- 
while mentioning at this point that the matching of these two solutions at r = 1 arises 
naturally in a replica saddle-point calculation of (^tr log (z* — <ft)(z — </>)) ~ F(0; r) 
along the lines of 0, as we now describe briefly. Before the replica limit n — > is 
taken, one deals with (det n (z* — <f>*)(z — <j))). This expression may be considered as 
the averaged fermion determinant in the fluctuating gluon field <fi when there are n 
fermion flavors degenerate in mass. The role of fermion mass is played by r = \z\. 
For r < 1 there are two large N saddle points, call them A and B, that govern the 
integral over certain replica matrix fields in the expression for the averaged deter- 
minant. These saddle points are obviously in 1 — 1 correspondence with the two 
solutions of ( |3.25D at rj = that we mentioned before. One of these saddle points, 
say A, dominates the other (i.e.; its integrand is larger.) However, when r > 1, the 
previously dominating saddle-point A disappears, and the determinant is given en- 
tirely by saddle-point B. The two expressions for this averaged determinant as well as 
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their first derivatives with respect to r match continuously at r = 1. In other words, 
the fermion "free energy" log (det n (z* — ft)(z — 0)) undergoes a second order phase 
transition at r = 1 as the fermion mass is varied. 
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4 Boundaries of the Eigenvalue Distribution 



In some applications of random non-hermitean matrices, one is primarily interested 
in the domains over which the density of eigenvalues is non-zero. In this section, we 
will develop a formalism for determining the boundary of these domains. 

First, let us get oriented. We have already remarked that the deterministic piece 
in H has N eigenvalues +r = \z\ and N eigenvalues — r. For small r = \z\ we would 
expect the random piece in H to smear these two adjacent delta function spikes in 
uj(fi; z, z*) into one continuous lump. In other words, the density of eigenvalues uj as 
a function of p is non-zero along a single segment. For large r, the deterministic part 
in H will dominate, and we expect uasa function of \x to separate into two lumps 
centered at ±r. On the other hand, in the original non-hermitean matrix problem, 
as \z\ increases, we expect to arrive at the boundary of the domain outside of which 
p(x, y) vanishes. 

It is thus natural and intuitive to expect that the criterion for determining the 
boundary in question is to locate the points in the complex z plane at which u(p; z, z*) 
splits into two lumps. Note that since u is an even function of p, as mentioned earlier, 
the split first occurs at \i — 0. We can readily formalize this expectation. Consider 
the positive semi-definite hermitean matrix 

h 2 = (z* - <f>i)( z - <f>) (4.1) 

whose eigenvalues we denote by hf. Then from fl3.7|) we have 

1 N 

*, **) = - ^) + ^ + I*)]) ■ (4-2) 

Thus, if w(0; z, z*) vanishes, that implies h and hence z — <fi (considered averaged over 
the randomness of course) do not have any zero eigenvalues, or equivalently, does 
not have any eigenvalue equal to z. This proves our expectation. 

How do we determine the boundary of the domains over which p(x, y), the density 
of eigenvalues of the non-hermitean matrices 0, is non-zero? To proceed, we will now 
broaden our horizons and discuss the problem for an entire class of which the Ginibre 
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problem that concluded the previous section is an example, namely the class in which 
the Green's function Q{r\\ {r}) of the associated hermitean problem is determined by 
the solution of a polynomial equation F(Q) — 0. In general, as we will see, there can 
be a whole set of parameters {r}; in Ginibre's problem there is only one real variable 
r. 

Let us decompose 

Q = u + iv (4.3) 

into its real and imaginary parts and imagine moving along just below the real axis 
(that is, we will take 77 = \x — ie as in the previous section.) Note that for rj real, the 
coefficients in the polynomial F(Q) are real as well. There will be segments over which 
v is non-zero; these segments are where the density of eigenvalues of the hermitean 
matrix H is non-zero. (Up to an irrelevant factor, v > is just a convenient notation 
for what we called u.) Let us now imagine moving along such a segment. We would 
like to determine the value of fi (call it fi c ) at which v vanishes, in other words, the 
endpoints of the segment over which v is non-zero. 

Let us decompose the equation F(Q) = into its real and imaginary parts and 
expand for v small: 

F(u) + h(u)v 2 + 0(v A ) = and 



v{{F'{u) + 0{v 2 )) = 0. (4.4) 

Here h(u) is a function whose detailed form does not concern us, and F' denotes the 
first derivative of F with respect to its argument. All these functions F, F', and h 
are of course also functions of fi (and of {r} which we think of as fixed.) 

At = [i c , v vanishes by definition and hence according to the first equation in 



(|4.4j) F(u) also vanishes, generically as F(u) ~ (fi — /i c ). Thus, we conclude that 
v ~ \fji — fjL c \^. (This type of argument is of course reminscent of the arguments 
used in Landau in his theory of phase transition. As a matter of fact, it explains 
the robustness of the square root singularity in u(/jl) near the edges of the eigenvalue 



distribution for many ensembles of hermitean matrices.) 

For [i away from /j c , v is non zero, and thus we may take out the factor of v from 
the second equation in ( |4.4j) and write F'{u) + 0{v 2 ) = 0. Letting fi go to /i c , we 
argue by continuity that F'{u) also vanishes at r\ c . 

Should the reader feel somewhat uneasy over the continuity argument, we can 
obtain the same conclusion by differentiating the second equation in ( |4.4| ) with respect 
to /i. Using v ~ \fi — [i c \? we have that as \i approaches /i c , ^ tends to infinity while 

tends to zero. We reach the same conclusion above. 

Ct/i 

An alternative argument is to write the polynomial F(Q) = Y\i(Q — Qi) in terms 
of its roots Qi. Since all coefficients in F(Q) are real (we assume that rj = fi — ie), 
the Qi are real or appear as complex-conjugate pairs. In particular, the "physical" 
root Q(n; {r}) is generically complex, and thus appears with its complex-conjugate 
partner. As /i — > fi c (from within the segment along which v > 0) v — > and thus 
the two roots Q(fi; {r}) and Q*(fi; {r}) collide, that is to say, become equal (and real, 
so that v (/i c ) = 0). Thus, F(Q) has generically a double zero at /x = /i c , and so both 
F and its first derivative F' vanish at that point. 

To summarize, to determine the boundary or edge of the density of eigenvalues of 
the hermitean random matrix H we solve the two equations F(u) = and F'(u) = 
to obtain function of the parameters {r}. In general, there will a whole set of 

This method gives us more than what we need to know in the context of the 
non-hermitean problem. To determine the boundary of the density of eigenvalues of 
the non-hermitean random matrix 0, we are interested in the point, as we vary {r}, 
when one of the /i c appears at the origin. Thus, we finally obtain the equations to 
determine the boundary in question, namely 



F(u) 



!a* c =o 



and 



F'(u) 



Mc=0 



(4.5) 
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To see how simply this procedure works apply it to Ginibre's problem in which 
F(u) = u 3 — 2r]u 2 + (1 + rf — r 2 )u — rj. We obtain almost instantly r = 1. Indeed, 
it is known from Ginibre's work, as well as from our derivation in Section 3, that the 
eigenvalues are distributed over a disk of radius unity. 

In all the examples discussed here and below, F(Q) is a polynomial of odd degree 
in Q. In Ginibre's problem F is a cubic, while in the examples below, F is a quintic. 
Furthermore, from (|3.15|) and the chiral structure of H, we know that Re Q{rf) is an 
odd function of rj. Thus, at 77 = \i — the real part of Q vanishes. In general, the 
imaginary part does not vanish of course, but precisely when the density of eigenvalues 
splits into two lumps, the imaginary part of Q also vanishes at \i — 0. It is clear 
from ( |3.15|) and the symmetry of u;(/x) that in the cut complex plane, away from the 
cut, Q{rj) is an odd function of its argument. Thus, if we write F(Q) = J2k c k{v)G k , 
the coefficient functions Cj(r]) must be even in 77 for odd j, and odd for even j, namely, 
obey 

Cj (- V ) = (-iy +1 Cj (r)) . (4.6) 

One can check by inspection of ( |3.25| ) that this holds in particular for Ginibre's 
polynomial. 

This suggests a computationally simpler procedure. Given the polynomial equa- 
tion F(Q) = 0, first set 77 to zero. Now from ( |4.6|) co(r? = 0) = 0, and so we should 
divide F(Q) by Q and set Q to zero (because we wish to identify the points {r} where 
\i c = 0.) At the risk of being repetitive but for the sake of clarity, let us apply this 
to Ginibre's polynomial Q 3 — 2i]Q 2 + (1 + rj 2 — r 2 )Q — r] — 0. After the first step, we 
have Q 3 + (1 — r 2 )Q = 0. After the second step, we have r = 1 as desired. 

The advantage of this procedure is that in the computation of Q from the equations 
analogous to ( [3.21|) and ( |3.23| ), we can from the beginning set 77 to zero and then 
compute various quantities to either zeroth order or first order in Q. This results in 
enormous simplifications. 
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5 Further Examples of the Application of the Method 
of Hermit izat ion 



The development of the method of hermitization in the previous two sections was 
formal and generic. In the last part of Section 3 we demonstrated the method of 
hermitization by applying it to Ginibre's problem. In this Section we provide further 
examples to the application of this method. 

5.1 Deterministic plus random, both non- hermit ean 

Consider a situation in which the fluctuating non-hermitean matrix <fi is shifted by 
a non-hermitean deterministic piece Hq. For simplicity, we take h to be diagonal, 
and also assume that <fi is drawn from the Gaussian distribution Q3.2C| ). Then our 
discussion of Ginibre's problem (Eq. (|3.21|) et. seq.) follows through with the obvious 



shifts z — > z — ho and z* — ► z* — h . In particular, ( J3.24 ) now reads 

/ Q -7] z - h 



=7 v 



R 2 -(V- Q) 2 



\z* -h\ g-v 



(5.7) 



where R 2 is the N x N diagonal matrix (z — hoy{z — h ). Setting rj = and tracing 
both sides of (|5.7Q , the "gap-equation" ( |3.25|) changes into 



Eq. ( |5.8| ) has a trivial solution Q = and a non-trivial solution which satisfies 
(l/iV)tr[l/(i? 2 — Q 2 )] = 1. Note that according to the analysis in Section 4, the 
boundary of the eigenvalue distribution occurs when the non-trivial solution of (|5.8|) 
vanishes as well, and is thus given by 

— tr — r^ = l- 5.9 



Extracting the lower left block of Q£ and tracing, we find from (|3.5| ) that 
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Substituting Q = into ( |5.10| ) we find p(x, y) = (1/N) tr S^^z — h ), which is simply 



the spectrum of eigenvalues of h , and it obviously vanishes outside the domain of the 
eigenvalue distribution of the random matrix. We thus associate the solution Q = 
(|5.8|) with the region of zero eigenvalue density outside the domain of eigenvalue 
distribution, as we already observed in our analysis of Ginibre's problem at the end 
of Section 3. The eigenvalue distribution is bounded inside the curve ( |5.9|) , and the 
density of eigenvalues is given by Q5.10|) with Q being the non-trivial solution of ( |5~8|) . 
It is instructive to work out some simple examples of this case. 

5.1.1 Circular deterministic spectrum 

Assume that the eigenvalues of h lie equally spaced along a circle of radius a. In the 
limit N —>■ oo they form a uniform continuum. Defining 

r dw 1 

I = J 2^1 W 2_r* + as-g* w+1 ' 
ar 

where w = exp i9 runs along the unit circle, we have for the non-trivial factor of (|5.8|) 



I=-ar. (5.12) 

Let us denote the roots of w 2 — [{r 2 + a 2 — Q 2 )/ar}w + 1 = by w\ and W2- Note that 
the roots of the quadratic equation satisfy W1W2 = 1. Thus, one of the roots, call it 
w in , always lies inside the unit circle, and the other, call it w out , outside. We obtain 



I = § 1 = 1 . Then (03) reduces to (a 2 + r 2 - Q 2 ) 2 - Aa 2 r 2 = 1. 

J 2m (w— wi)(w— wo) w ir ,—w ou t v l V v •'J 



2-ni (w-wi)(w— V12) w in —w out ' 

There are obviously two roots for Q 2 , but one of them is positive for all values of r, 
and thus cannot correspond to the eigenvalue distribution. The remaining (physical) 
root is 



g 2 = r 2 + a 2 - Vl + 4a 2 r 2 . (5.13) 
The boundary of the eigenvalue distribution occurs at Q = 0, namely at 

r 2 = a 2 ± 1 . (5.14) 



Thus, if a > 1, the original circle of eigenvalues swells into the annulus va 2 — T < 
r < Va 2 + 1, while for a < 1, the eigenvalues fill a disk of radius \/a 2 + 1. Finally, 
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from ( |5.10D we find the density inside the annulus (or the disk) as 



p(x,y) = - (1 ; a 1 . (5.15) 

HK ,yJ Ti V VI + 4a 2 r 2 J V ; 

Note that there are more eigenvalues near the outer edge. For the annulus case (a > 1) 
the densities at the outer and inner edges are given by 7rp(outer) = 1 — a 2 /(2a 2 + 1) 
and 7rp(inner) = 1 — a 2 /(2a 2 — 1). For the disk case (a < 1) the density at the origin 
is given by 7rp(0) = 1 — a 2 and is thus reduced from the density in Ginibre's case by 
a factor of 1 — a 2 . 

In this example the eigenvalues of the deterministic piece ho were located along 
a circle. It would be interesting to generilize this example to situations where the 
eigenvalues of ho are located along an arbitrary curve in the complex plane. 

5.1.2 Two point deterministic spectrum 

In our next example we consider an ho with eigenvalues ±e, equally degenerated. 
With no loss of generality we may always take e real and positive. The non-trivial 
factor of ( |5.8| ) is now 

+ ... ■ ^ = 2 - ( 5 - 16 ) 



\z-e\ 2 -G 2 \z + e\ 2 -G 2 
Setting Q = in (|5.16f) we conclude that the boundary is determined by 



y 4 + (2x 2 + 2e 2 - l)y 2 + x 4 - (2e 2 + l)x 2 + e 4 - e 2 = . (5.17) 

Clearly, the eigenvalue distribution is an even function in both x and y. For e = 
this is simply Ginibre's problem and the eigenvalues are distributed uniformly on a 
disk of unit radius centered around the origin. As we increase e from zero, the disk 
is distorted: it gets squeezed along the y axis and stretched along the x axis (in both 
directions). We find that for e < 1 the eigenvalues are still concentrated in a single 
blob centered at the origin and confined within x 2 < (2e 2 + 1 + y/1 + 8e 2 )/2. However, 
at e = 1, what used to be the disk is now completely pinched and separates into two 
identical lobes that touch only at the origin. As we increase e above 1, these two lobes 
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separate and their extent along the x axis is determined by (2e 2 + 1 — vl + 8e 2 )/2 < 
< (2e 2 + l + vT+8e 2 )/2. The expression for p(x, y), (Eq. ( |5.10|) ) is not particularly 
illuminating and we do not bother to write it here. More generally, if ho has K 
different eigenvalues e^, Ginibre's disk would get torn into K different blobs as we 
scale ho up, obviously. 

5.2 Two point deterministic spectrum plus hermitean ran- 



As a complication of the previous subsection, we will now discuss the problem in which 
4> = h + M consists of the sum of a deterministic piece h and an hermitean random 
matrix M which we will take from the Gaussian distribution P(M) = tr m2 ■ 

This case is more complicated than the previous cases because the hermiticity of M 
amounts to a constraint on the "gluon propagator". Our interest in this problem is 
partly motivated by the flux line pinning problem to be discussed in Section 6. In that 
problem the eigenvalues of ho trace out an ellipse. Here we content ourselves with the 
much simpler case of ho being a diagonal matrix, with half its eigenvalues equal to it 
and the other half equal to — ie (with e real and positive.) The qualitative features 
of the eigenvalue distribution is clear. For e = 0, the eigenvalues are restricted to the 
real axis and follow Wigner's semi-circle distribution. As e increases, the eigenvalues 
invade the complex plane, and eventually, as e increases past some critical value e c we 
expect the domain over which the density of eigenvalues is non-zero to break up into 
two "blobs" centered around ±ze. Note the obvious symmetries of the shape of the 
eigenvalue distribution: if (x, y) is a point on the boundary, so are (±x, ±y). Here we 
shall content ourselves with calculating the boundary of the eigenvalue distribution 
only. Due to the specific structure of ho we break M into (N/2) x (N/2) blocks 



dom 



(AC 

v C t B 




(5.18) 



in terms of which the Gaussian distribution reads 



1 -f tr (A 2 +B 2 )-N tr CtC* 

z 



(5.19) 
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(In general, of course, we can consider taking the Gaussian widths of A, B, and C 
to be different, but we will not entertain this possibility here.) With and M given 
above, the deterministic piece of the Hamiltonian (371) in terms of its {N/2) x (N/2) 
blocks is 

/ ie-z \ 



H = -g- 1 ( v = o) 



o 



— ie — z 











— ie — z 



and the random part is 



V 



) 16- Z 

( A C \ 

Ct 5 

A C 



(5.20) 



B 



(5.21) 



Speaking picturesquely, and adding to the "quark-gluon" language mentioned earlier, 
we can refer to the (block) four dimensional basis as "flavor" , with flavor assignments 



(u\ 

c 

d 

V 8 J 



(5.22) 



containing two "generations", (u,d) and (c, s). With this flavor assignment, we see 
from the interaction ip^Vifj that the gluons corresponding to A, B and C all act like 
weak interactions gauge bosons. A and B act within their appropriate generation, 
while C and communicate between the two generations.^ Since we are only inter- 
ested in determining the boundary of the eigenvalue distribution, from now on, we 

3 To see this more clearly, and also for computational simplicity, it is convenient to change the 
basis by interchanging rows as well as columns 2 and 3. In the new basis the A's in V will migrate 
to the upper left corner, the B's to the lower right corner, the C's to the upper right block, and the 
C^'s to the lower left block, such that in the new basis V will become simply M ®<j\. However, we 
do not use the transformed basis in the text explicitly, to avoid unnecessary repetitions of formulas. 
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set 77 = 0, so that Q^(t] = 0) is hermitean. From ( |3.22| ) and ( |5.19| ) we then find that 
the self-energy matrix has the form 

/ #33 + 944 5-31 + g 42 \ 



#33 + 944 5f 3 i + 5f 42 
£31 + 9*42 9n + 922 

#31 +#42 #11+322 / 



(5.23) 



where the non-zero blocks are proportional to Ijv with proportionality coefficients 
that are sums of 



9ab — 9 ba 



N 



tr(S) 



ab block j 



a, b = 1, • • • 4 



(5.24) 



where the trace is carried over the (N/2) x (N/2) blocks. By "index democracy" 
and the texture of H, we obviously have gn = (722 and #33 = #44. Then clearly 
Q = (#11 +5'33)/2. Substituting (|5.23| ) into ( [3.21D we obtain a system of coupled 
"gap equations" for the coefficients gat,. We are interested only in the g a b that appear 
explicitly in £.[] Introducing the quantity V = (g n — g^)/2 we obtain the equations 
for these elements as 



V 
Q 

913 
924 



G(\9i3\ 2 -\924\ 2 -^V) 

G[\gi3\ 2 + \924\ 2 -2(g 2 + V 2 )] 

(\9i3\ 2 - G 2 )(z -it- g* 13 - g* M ) 
(\924\ 2 - Q 2 )(z + ie- g* 13 - g* 4 ) 



(5.25) 



These are 6 real equations in 2 complex plus 2 real unknowns. Being interested only 
in determining the boundary of the eigenvalue distribution, we only need the non- 
trivial factor of the gap equations to first order in (?(0). From the first equation in 
(|5.25|) we find that T> = when Q = (for the trivial, as well as for the non-trivial 

4 Note a lso t hat the trace of the lower left block of is g^i + .942- Thus the g a b that do not 
appear in (5.23) arc also irrelevant for the purpose of calculating p(x,y) from (2.11). 
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factor of the gap equation), thus eliminating one real unknown and one real equation. 
Then, dividing the second equation in ( |5.25| ) by Q and setting Q = we obtain 



0* -He - #13 - # 24 )#i3 

{Z* - it- #13 - # 24 )#24 



(5.26) 



which constitute 5 real equations in 2 complex unknowns, thus producing a single 
real constraint, namely, the desired equation for the boundary. A direct elimination 
of one of the g's in (|5.26|) produces a cubic equation in the other #, which has to be 
solved explicitly. Then, the correct root of that cubic has to be fed into the constraint 
to obtain the boundary. We will content ourselves with studying perturbatively the 
two limits e << 1 and e >> 1. 

We first examine the small e limit. It is useful to introduce the combinations 



#13 + #24 = gV2 , #13 - #24 = 



and the parametrizations 



2^2 sin/3, ie = \V2 . 



(5.27) 



(5.28) 



Then the appropriate linear combinations of the equations Q5.26 ) become 



# 2 - 2# sin - \f + 1 = 
gf - 2f sin - Xg = 



l/| S 



\g\ 



i . 



(5.29) 



At e = we should recover the ordinary Wigner semi-circular behavior associated 
with flSTTSD 

(o) (o) z* - Vz* 2 - 8 
013 = 024 = ^ (5-30) 

namely, 



i(o) = z . _ V ^2 _ 8 = ^ ^ = ^ 



2V2 
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Thus, in a systematic expansion of ( |5.29| ) in powers of A around A = 0, / is always 
a small quantity, which is the reason for introducing the combinations (|5.27f) . The 
eigenvalues lie on the real axis in the segment \x\ < 2\[2 which is the (degenerate) 
boundary of the eigenvalue distribution. In particular, along this "boundary", the 
angle (3 = (3^ (0 < (3^ < 2n) is real, of course, in accordance with the third 
equation in (|5.29| ). Let us now solve the first two equations in (|5.29 ) to leading order 
in A. We find 

9 = ^ f-1 + £*L) , / = Ae 2 ^ . (5.32) 
\ 2 cos [3 J 

Substituting these quantities into the third equation in (|5.29 ) we find that to this 
order 

(3 = p(°) - zA 2 , (5.33) 
and thus, we find from (|5.28|) that to this order the boundary is 



2 



x 2 + ^- = 8, (5.34) 

which shows that the original segment containing the eigenvalues of M swells into an 
ellipse, with a semi-minor axis equal to 2y / 2|A| 2 = \^2e 2 . 

We now turn to the other extreme limit, e >> 1. Here we expect of course that 
the eigenvalues concentrate in two semicircular distributions around z = ±ze. For 
concreteness, let us concentrate on the vicinity of z = —it. Shifting z* = ie + (, we 
find that the last two equations in (|5.26|) become 



924 + — = C - 013 > ( 2 «e + —)gi3 = 1 . (5.35) 

924 924 

We now expand (|5.35| ) in powers of 1/e, for |C|/e << 1- It is clear that here #13 is the 
small quantity. To leading order in 1/e we find the semi-circular behavior 

(0) C - VC 2 - 4 (!) 1 

9\l = g , 9 { J = ^ • (5.36) 

Thus, to this order in 1/e, the eigenvalues follow semi-circular distribution along 
the segments \x\ < 2,y — ±e, which are also the (degenerate) boundaries of the 
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eigenvalue distribution. Note that the width squared of this semi-circle is half that 
of the semi-circle in ( p. 31] ), because the eigenvalue distribution is now split into two 
identical blobs. It is natural to parametrize 

C = 2sin 7 , (5.37) 

in terms of which gf} = —ie n . Since ( is real along the "boundaries" to leading 
order in 1/e, 7 = 7^ (0 < 7^ < 2ir) is real as well, in accordance with the first 
equation in ( |5.26| ). Expanding g 2 4 to the next order in 1/e we find 

Substituting these expressions for g 2 4 and #i 3 into the first equation in (|5.26|) we find 
to this accuracy that on the boundary 

7 = 7 (0) - , - (0) , (5-39) 
4e cos 7 W 



and thus from ( 5.37p the (lower half of the) boundary is given parametrically by 

{ r=-' n + f « 5 - 4o) 

We see that the effect of the first non-trivial correction is simply to pull the two 
semi-circle distributions towards the real axis in an amount proportional to 1/e, with 
no distortions. We shall present a simple explanation of this phenomenon in Section 
6. 

We close this subsection by determining the critical value e c at which the connected 
eigenvalue distribution corresponding to small values of e breaks into two separate 
lumps. From the symmetries of the eigenvalue distribution it is clear that when 
this happens, the boundary passes throught the origin (it then has the shape of the 
figure 8). Thus, substituting z = (3 = in ( |5.29| ) we find a non-trivial solution 
g = 0, / = 1/A, and thus e c = y2. 
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5.3 QCD at finite chemical potential 



As was mentioned in the Introduction, StephanovpJ has recently studied the chiral 
phase transition in QCD at finite chemical potential using non-hermitean random 
matrices. The Euclidean Dirac operator at chemical potential \i is non-hermitean. 
The matrix representation of this operator (for massless quarks) is 

/ C + in 



(5.41) 







The chiral components of the Dirac operator are represented here by the N x N 
complex matrices C and C\ The fluctuations of the gluonic background are modeled 
by taking C to be Gaussian random, namely, with probability distribution P(C) = 
exp(—NtrC^C)/Z. The Green's function and density of eigenvalues of ( 5.41Q were 
calculated in || using the replica method, and were rederived in diagrammatically. 
Thus, as a complementary calculation, in the following we shall apply the methods of 
Section 4 to calculate the boundary of the eigenvalue distribution for this ensemble, 
without calculating the density itself. Recall that to the sole purpose of determining 
the boundary, we only need the the non-trivial factor of the gap equation to first 
order in Q at rj = 0. Thus, as in previous sections, we set rj = from now on. With <p 
given by ( |5.41| ), the deterministic piece of the Hamiltonian ( |3.1|) in terms of its N x iV 
blocks is 



h = -g^( v = o) 



( -z ifi \ 

ifi —z 

-z* -ifi 

i -ifi -z* J 



(5.42) 
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and the random part is 



V 



( c \ 
o 
C 



(5.43) 



\ ct j 

Then from ( |3.22| ) we findF] that the self-energy matrix has the form 

I #44 5f 42 \ 

#33 #31 
#42 #22 
V #31 #H J 

where the non-zero blocks are proportional to l^v with proportionality coefficients 



(5.44) 



9ab — 9ba 



block 



a,6 = 1, • • -4. 



(5.45) 



By "index democracy" all the diagonal coefficients gkk are equal, and are therefore 
equal to Q. Similarly, g 42 = g 31 = G (they are in the same block of E.) Substituting 
(|5.44|) into (|3.21|) we obtain a system of coupled "gap equations" for the coefficients 
# a fe. We are interested only in the g a b that appear explicitly in E.Q Among the gap 
equations there is subset of one real equation and one complex equation that involve 
only the real unknown Q and the complex unknown G. In accordance with our 
discussion in Section 4, the real equation contains Q as an overall factor. Following 
the steps indicated in Section 4 we divide this factor out and set Q = 0, thereby 
obtainig one equation for G and a constraint. Eliminating G, we find the constraint 

/i 2 (4/i 2 + x 2 ) [Gu 2 - y 2 fx 2 + 4/iV + (1 + 4/i 2 - 8//)# 2 + 4/. V - 1)] = . (5.46) 



5 See footnote 3. Here of course, we do not have the blocks A and B. 

6 Note also that the trace of the lower left block of G„ is (331 + <?42)/2 = G. Thus the g a b that do 



not appear in (5.44) arc also irrelevant for the purpose of calculating p(x,y) from (2.11) 



31 



Solving (|5.46|) for x 2 we obtain 
1 



x 



(/i 2 -y 2 ) 2 
in accordance with [Efl. 



(5.47) 
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6 Energy Level Dynamics of Non-Hermitean Hamil- 
tonians 



In much of this paper, we have discussed non-hermitean random matrices without 
any spatial structure. Physical problems typically involve spatial structures and thus 
random matrices with many zero entries, with non-zero entries representing nearest 
neighbor hopping for instance. For such problems, even in the hermitean case, there is 
typically no simple analytic approach available and these problems has to be studied 
on a case by case basis. Here we will use the problem of driven vortex lines as an 
illustrative context within which to make some general qualitative comments. In the 
simplest model for this problem, the Hamiltonian H is the sum of a deterministic 
but non-hermitean hopping matrix Hq and a random but hermitean matrix repre- 
senting pinning by impurities V. We have in the simplest one-dimensional case with 
iV sites 

H Wj = e h 5 i+1J + e~ h 5 iJ+1 , i,j = l,---N (6.1) 

describing a particle which tends to hop in one direction more than in the opposite 
direction. The eigenvalues of H (assuming periodic boundary conditions) can be 
easily worked out by standard methods to be 

27TTi 

E n = 2 cos (— -ih) {n = 0, 1, • • • , N - 1) , (6.2) 

thus tracing out an ellipse. (The limits of h — > and h — > oo in which the ellipse 
collapses to a line segment on the real axis and expands to a circle are both easily 
understood.) The precise form of V reflects what one believes to be an appropriate 
representation of the microscopic physics, of course. Some choices that come to mind 
are a random nearest hopping Hamiltonian = Vi + v* <5ij+i, or a diagonal 

Hamiltonian describing site impurities Vy = V{ 5ij, or a purely random hermitean 
matrix V. This last case may be most amenable to analytic study, but represents the 
somewhat unphysical situation of being able to hop from any site to any other site 
(the "infinite dimension" limit). 



33 



Let us begin our discussion of random non-hermitean hopping by making some 
simple observation based on symmetry. For H real, the non-real eigenvalues come 
in complex conjugate pairs. This is true in particular for H in ( |6.1| ), namely E n 
and E^-n in (|6.2|) form such a pair. Furthermore, if ip is an eigenvector (say, on the 
right) of Hq, corresponding to an eigenvalue E, ip* is an eigenvector (on the right) 
corresponding to an eigenvalue E*. For the specific Ho in ( |6.1| ), if N is even, then the 
eigenvalues also come in opposite pairs: E n = —E n+ N/2- This follows for a large class 
of hopping Hamiltonians with "bipartite symmetry" under an operation in which we 
multiply the basis vectors associated with the odd sites by (—1). 

As we crank up the hermitean V we expect the eigenvalues tracing out the ellipse 
to slowly migrate to the real axis, and becoming real in the infinite V limit. In Figure 
(1) we present numerical results for the spectrum of H + V (with periodic boundary 
conditions), where V describes random site energies. These site energies are taken 
from a flat distribution, symmetric about the origin, and their range is comparable 
with the size of the ellipse. We do not bother to display here how the distortion of 
the ellipse varies with the range of the random site energies. 

To the initial surprise of the authors, already for a finite amount of randomness, 
there are pairs of eigenvalues which have "snapped" onto the real axis. In fact, this 
"snapping phenomenon" is easy to understand qualitatively. 

Repulsion of nearby energy eigenvalues is of course one of the fundamental features 
of quantum mechanics, and is explained in any elementary texts as a property of 2 x 2 
hermitean matrices. If we have a pair of eigenvalues ±e, and turn on an off-diagonal 
perturbation v between them, the eigenvalues become ±y/e 2 + v 2 , thus moving apart. 
The behavior of random hermitean matrices can be understood qualitatively in this 
way, and is formalized in the concept of the Dyson gas with its molecular repulsion. 

Now if we have a pair of complex eigenvalues separated in the imaginary direction, 
that is, if we have a pair of eigenvalues ±ie, we see that under the (hermitean) 
perturbation v they become ±i\/e 2 — v 2 . Our standard intuition fails. Two nearby 
imaginary eigenvalues attract! Furthermore, they do not migrate gradually to the 
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0.75 - 
0.5- 

0.25 - 



-0.25 - 
-0.5 - 
-0.75 - 



Figure 1: The spectrum of ( |6.1| ) with random site energies (with 100 sites). Here 
h = 0.5 and the random site energies \vi\ < 1.5 are drawn from a flat distribution, 
symmetric about the origin. There are 12 real eigenvalues in all. 
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real axis. Rather, as soon as v becomes larger than e in magnitude, the eigenvalues 
snap to the real axis. 

(We note that in the case of an anti-hermitean perturbation (that is, for v imag- 
inary) the opposite behavior holds: two nearby eigenvalues separated along the real 
direction attract, and two nearby eigenvalues separated along the imaginary direction 
repel.) 

As in the case of hermitean matrices, we can reaffirm the mutual effects of nearby 
energy levels we observed in the simple 2x2 matrix example by doing a second order 
perturbation calculation. For simplicity, we assume that Hq is non-degenerate. Let 
us then add a perturbation V to Hq. We recall that a non- hermitean Hamiltonian H 
with eigenvalues {E^} has a set of eigenvectors on the right H \ip^) = E^\ip^), 
as well as a set of eigenvectors on the left (4>^\H = (<fi$\E%\ that form a complete 



bi-orthonormal system[23|. With these facts in mind, perturbation theory for non- 
hermitean Hamiltonians appears formally the same as in the purely hermitean case 
(as long as one remembers to make the distinction between eigenvectors on the right 
and eigenvectors on the left.) In particular, the correction of E^> to first order in V 
is simply 5^' E n = ((f>^\V\ip^'}, and the correction to second order in V is 

e (2) F _ v mvi^i^ivm ( , 

For a diagonal Hq, the bi-orthogonal sets coincide ipffl = 4>^\ forming a single or- 
dinary orthogonal basis, and (|6.3|) assumes the same form as in perturbation theory 
of hermitean operators. In this case we see immediately that if V is hermitean, the 
original eigenvalues repel in the real direction and attract in the imaginary direction, 
whereas if V is anti-hermitean, they do the opposite, as expected. 

To get a qualitative feel, let us analyze the effect of a pair of eigenvalues on 
another nearby pair. We concentrate on the 4x4 piece of the Hamiltonian associated 
with these eigenvalues, ignoring all the other eigenvalues and take H = Hq + V with 
Hq = diag(ie, — ie, c + ia,c — ia) ; e, a > 0. For simplicity, we take the perturbation V 
hermitean with zero diagonal entries, and all other entries real and equal to v (this 
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should be enough to capture the qualitative behavior in a typical case.) Assuming 
that e << |c| ~ a, the pair of eigenvalues of H corresponding to ±ie would snap first, 
and we want to study the effect of the pair at c±ia on it. Cranking up ( |6.3| ) we thus 
find that to second order in perturbation theory, the first two eigenvalues are 



1 1 

+ 



c + i(a — e) c — i(a + e) 



which to first order in e/c 2 reads 



/ v 2 9 c 2 — a 2 \ lev 2 _ / 9 9 . ,\ 

± !e i 1 " 2? ~ 2v wr^r )-^ + (" e /a ) ' (6 - 4) 

We readily indentify the piece ±ie (l — in (|6.4|) as the first two terms in the expan- 
sion of iix/e 2 " — v 2 we mentioned earlier, namely, the exact value of these eigenvalues 
when c ± ia — > oo. Thus, under the influence of V, ±ze indeed attract each other. 
Note that to the order to which ( |6.4j ) holds, the pair c±ia slows or accelerates the 
snapping of the first pair onto the real axis depending on whether c 2 < a 2 or c 2 > a 2 . 
The first pair is also repelled by the pair c±ia along the real axis. If c > it moves 
to the left, and if c < it moves to the right. 

The elementary discussion just concluded explains qualitatively the distortion of 
the ellipse corresponding to Hq that we see in Figure (1). The eigenvalues near the 
tips of the ellipse are the first to snap to the real axis, for the obvious reason that 
they are closer to the real axis. They are then repelled away from the ellipse by 
the eigenvalues that remain off the real axis. Once on the real axis, they repel the 
eigenvalues that snapped earlier further along the real axis. 

For random hermitean hopping Hamiltonians, the lore is that the states in the 
middle of the band (or possibly only one state) are extended, while the states on 
the two "wings" of the band are localized. For the one- dimensional case, this can be 



shown quite explictly by analytical methods ||25||, and for the two-dimensional case, 
this can be studied by renormalization group arguments and verified numerically. In 
the present one dimensional context, we might thus expect (or conjecture) that the 
eigenvalues that have snapped onto the real axis correspond to localized states, while 
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Figure 2: The eigenvector ipi associated with one of the eigenvalues (2.9) on Fig. (1) 
that snapped onto the real axis. The corresponding normalized participation ratio is 
0.04. 



the eigenvalues that remain complex correspond to extended states 0. In Figures 
(2) and (3) we present numerical studies of the participation ratio that support this 
conjecture or expectationF]. The normalized participation ratio for a state ipi is defined 

as 



p 



(6.5) 



and thus should tend to 1 for extended states and to for localized states. 

It is worth emphasizing that the method of hermitization given in Section 3 is 
independent of the spatial character of the random non-hermitean Hamiltonian. Thus, 
in principle, it may be used to reduce the problem of random non-hermitean hopping 



7 Out of the 12 real eigenvalues we found numerically in this case, we chose to depict the most 
localized state. Some of the other states corresponding to real eigenvalues may break into 2 or 3 
separated lumps, but their participation ratios remain of the same order of magnitude as of the state 
presented on Fig. (2). 
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Figure 3: The eigenvector associated with one of the complex eigenvalues (0.2 + 
0.8 i) on Fig. (1) . The corresponding normalised participation ratio is 0.8. 
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to the auxilary hermitean problem. The deterministic part of the hermitean problem 
can be interpreted as a nearest neighbor hopping of a particle with a binary internal 
state (call it an "up quark" or a "down quark"), such that whenever the particle hops 
it flips its internal state. In any case, it still bears the same spatial structure as in 
the original non-hermitean problem, thus retaining all the difficulties associated with 
random hermitean matrices with spatial structure. 
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Appendix : A 2 x 2 example 



In this Appendix we will perform an exceedingly elementary but instructive exer- 
cise. Consider the 2x2 non-hermitean matrix 

C + ip 
C* + ip 

where C is a Gaussian random complex variable and p is a real parameter. This is 
thus the N — 1 limit of the problem considered by Stephanov||, which we discussed 
in sub-section 5.2. The eigenvalues of <p are A = ±[(C + ip)(C* + ip)] l l 2 . Writing 
C = re %e we have by definition the density of eigenvalues 

p{x, y) = J J dr 2 ^e~ r2 5(x - Re A) 5{y - Im A) (A.l) 

The evaluation of this integral is elementary and gives 

p(x, y) = - , (x2 + y2) e^ 2 -y 2 ^ 9(p 2 - y 2 ) . (A.2) 



^ v /(/i 2 + x 2 )(/i 2 - ? / 2 ) 



As is indicated by the 9 function, the density of eigenvalues is non-zero over an 
infinitely long strip with width 2p. It decreases exponentially as x —>■ ±oo and blows 
up as an inverse square root as y — > ±/i. Notice that if we impose the condition that 
H be normal (but non-hermitean), that is [0, (jy] = 0, C is required to be real, and 
the density of eigenvalues collapses to 

p(x, y) = l -P{x) [5(y - p) + 8{y + p)] (A.3) 

where P{C) denotes the (normalized) probability of distribution of C, whatever one 
might choose. The infinitely long strip collapses into two lines. 
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